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Abstract 

We argue that many aspects of improper ferroelectric activity in manganites with the Pbnm and 
P2\nm orthorhombic structure can be rationalized by considering the limit of infinite intra-atomic 
splitting between the majority- and minority-spin states (or the double exchange limit), which 
reduces the problem to the analysis of a spinless double exchange (DE) Hamiltonian. We apply this 
strategy to the low-energy model, derived from the first-principles electronic structure calculations, 
and combine it with the Berry-phase theory of electric polarization. We start with the simplest two- 
orbital model, describing the behavior of the e g bands, and apply it to the E-type antiferromagnetic 
(AFM) phase, which in the DE limit effectively breaks up into one-dimensional zigzag chains. We 
derive an analytical expression for the electronic polarization (P cl ) and explain how it depends on 
the orbital ordering and the energy splitting A between e g states. Then, we evaluate parameters of 
this model for the series of manganites. For these purposes we start from a more general five-orbital 
model for all Mn 3d bands and construct a new downfolded model for the e g bands. From the 
analysis of these parameters, we conclude that the behavior of P cl in realistic manganites always 
corresponds to the limit of large A. This property holds for all considered compounds even in the 
local-density approximation, which typically underestimates A. We further utilize this property 
in order to derive an analytical expression for P el in a general two-fold periodic magnetic texture, 
based on the five-orbital model and the perturbation-theory expansion for the Wannier functions 
in the first order of 1/A. This expression explains the functional dependence of P el on the relative 
directions of spins. Furthermore, it suggests that P cl is related to the asymmetry of the transfer 
integrals, which should simultaneously have symmetric and antisymmetric components. Finally, 
we explain how the polarization can be switched between orthorhombic directions a and c by 
inverting the zigzag AFM texture in every second ab plane. We argue that this property is generic 
and can be realized even in the twofold periodic texture. 
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I. INTRODUCTION 



The multiferroic materials (or multiferroics) , where ferroelectricity coexists with some 
long-range magnetic order, have attracted a great deal of attention.- A very special class of 
multiferroics is improper ferroelectrics. In the latter case, the ferroelectric (FE) polarization 
not only coexists, but can be induced by the magnetic order. The improper ferroelectrics 
are expected to display a strong magneto-electric coupling, which is extremely important for 
practical applications. For instance, because of such coupling, the FE polarization can be 
efficiently controlled by the magnetic field, while the magnetization can be controlled by the 
electric field. From a technological point of view, the ultimate goal is to find materials with 
the large FE polarization, which would be coupled to the magnetic texture at maximally 
possible temperature (meaning that the magnetic transition temperature should be also 
high) . 

Manganites, crystalizing in the orthorhombic Pbnm and P2inm structure, are regarded as 
one of the key multiferroic materials. Despite low magnetic transition temperature (typically, 
less than 40 K) and modest values of the FE polarization (less than 1 /zC/cm 2 ), which have 
been achieved so far,- they have all essential ingredients to be called improper ferroelectrics. 
Namely, the appearance of ferroelectricity coincides with some complex magnetic ordering. 
Moreover, the possibility of switching the electric polarization by the magnetic fields has 
been directly demonstrated experimentally- Therefore, these materials are fundamentally 
important and are typically used as a playground for testing various theories and models of 
multiferroicity. 

Nevertheless, the theoretical understanding of improper ferroelectricity in these com- 
pounds is still rather controversial and there is no unique view on the origin of this effect. 
First, all multiferroic manganites are rather artificially divided in two groups: 

(i) the systems with the twofold periodic E'-type antiferromagnetic (AFM) texture (such 
as HoMnOa and YMnOa), where the FE activity is attributed to the nonrelativistic 
exchange striction,-^ and 

(ii) the rest of the systems, with more general magnetic periodicity, where the FE activity 
is believed to be due to the relativistic spin-orbit (SO) interaction and the magnetic 
texture itself is ascribed to the spin spiral.- The typical example of such systems is 
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TbMnOs, which has nearly fourfold periodic magnetic texture. 

This point was rationalized in the previous publications of one of the authors (Ref. |7| and|8|), 
where it was argued that there is no conceptual difference between twofold periodic and other 
multiferroic manganiets. The relativistic SO interaction plays an equally important role in 
both it deforms the E-type AFM state in the direction of the spin spiral, it will 

also deform the spin spiral and form a more general spatially inhomogeneous magnetic state. 
Thus, the ground state of multiferroic manganites will be neither the collinear Estate nor 
the homogeneous spin spiral. The relativistic SO interaction is essential for producing this 
inhomogeneity. However, the FE polarization itself is a nonrelativistic quantity in the sense 
that, for a given inhomogeneous distribution of spins, the appearance of the FE polarization 
can be described by nonrelativistic theories. 

Another group of controversies is related to the question: How to calculate the polariza- 
tion and what is the main contribution to it? Most of model calculations rely on the purely 
ionic picture, where the noncentrosymmetric distribution of spins gives rise to noncentrosym- 
metric atomic displacements. Then, the polarization is evaluated in the framework of the 
point charge model.-^ On the other hand, all modern first-principles calculations of the FE 
polarization are based on the Berry-phase theory.—^ Besides ionic polarization, the Berry- 
phase theory prescribes the existence of an electronic term. The latter can be expressed 
through the Wannier functions and is reduced to the ionic polarization only if the Wannier 
functions are fully localized at the atomic sites. In this sense, the deviation from the ionic 
picture is a measure of itineracy of the system. Moreover, unlike the ionic contribution, the 
electronic polarization can be finite even in the centrosymmetric crystal structure, provided 
that the inversion symmetry is broken by a magnetic order. Thus, the Berry-phase theory 
excellently suits for improper ferroelectrics. The first-principles calculations show that the 
electronic polarization can be as large as or even exceed the ionic contribution.- Neverthe- 
less, the physical meaning of this effect is still rather obscure and the electronic polarization 
is largely ignored in model calculations of multiferroic manganites. 

The purpose of this work is to make a bridge between first-principle electronic structure 
calculations and models of the FE polarization. Our main message is that the electronic po- 
larization is important and cannot be ignored. In the model calculations, it can be described 
by some "superexchange type" theories, similar to interatomic magnetic interactions.— 1 ^ On 
the other hand, in the first-principles calculations, one should pay a special attention to the 



relative direction of the electronic and ionic polarization: because of additional approxima- 
tions, results of theoretical structural optimization do not necessarily guarantee the correct 
answer to this question. 

Our analysis will be based on results of two previous works (Refs. Q and fj), where 



(i) A realistic low-energy model for the Mn 3d bands of manganites was constructed 
on the basis of first-principles electronic structure calculations in the local-density 
approximation (LDA); 

(ii) This model was applied for the search of the magnetic ground state of orthorhombic 
manganites; 

(iii) The model calculations were supplemented with the Berry-phase theory for the analysis 
of the FE polarization and its dependence on the form of the magnetic ground state. 

In this work we will further rationalize the story. First, we will show that the behavior of 
the FE polarization can be well described in the framework of the double exchange (DE) 
theory— The definition of the DE Hamiltonian will be given in Sec. [Ill Particularly, we 
will show that with the proper definition of the DE model, which should include effects of 
orbital polarization of Coulombic origin, one can reproduce, even quantitatively, the values 
of FE polarization obtained in a more general mean-field Hartree-Fock (HF) calculations 
for the low-energy model. Then, we will introduce an analytically solvable model for the e g 
electrons in the single zigzag chain (Sec. IIII A[) and argue that, besides double exchange, the 
behavior of electronic polarization in realistic manganites always corresponds to the limit of 
large intra-atomic energy splitting A between e g states (Sec. IIII Bl) . It will allow us to further 
generalize our story and derive an analytical expression for the electronic polarization in an 
arbitrary twofold periodic magnetic texture, based on the perturbation theory expansion for 
the Wannier functions in the first order of 1/A (Sec. IIII C]) . The idea itself has some similari- 
ties with the superexchange theory of interatomic magnetic interactions.— 1 ^ This analytical 
expression nicely explains the behavior of electronic polarization in the low-energy model as 
well as in the more general first-principles calculations. It also provides a good quantitative 
estimate for the polarization. In Sec. IIII Dl we will present a critical analysis of relative 
directions of electronic and ionic polarizations in the experimental and theoretically opti- 
mized P2inm structures of YMnC>3. Then, in Sec. IIII we will explain how the electronic 



polarization can be manipulated by changing the magnetic texture. Finally, in Sec. IIV} we 
draw our conclusions. 

II. BASIC IDEA AND APPROXIMATIONS 

The starting point of our work is that the main electronic and magnetic properties of 
multiferroic manganites can be described reasonably well by the one-electron Hamiltonian: 

H^ = k 3 + VA V (1) 

which is constructed in the basis of Wannier orbitals for the Mn 3d bands. In this notations, 
the matrix t^j has site-diagonal (i = j) and off-diagonal (i 7^ j) elements: the former 
describes the crystal-field effects, while the latter stands for transfer integrals. We do not 
consider explicitly the relativistic spin-orbit (SO) interaction. More specifically, it is assumed 
that the SO interaction is important for specifying the directions of spins in some noncollinear 
magnetic texture. However, it is unimportant for calculations of the FE polarization itself, 
provided that the directions of spins are known and the corresponding magnetic texture can 
be described by appropriate rotations of the mean-field potentials Vj, which will be specified 
below. Therefore, the matrix Uj does not depend on the spin-indices, s(s')= t or I, and 
can be presented in the form iij = ||£^ m '5 ss /||. In the more general five-orbital model, that 
we consider, the indices m and m' have the following order: m(m')= xy, yz, 3z 2 —r 2 , zx, or 
x 2 —y 2 . In the two-orbital model, constructed only for the e g bands, the indices m and m' 
run over 3z 2 — r 2 and x 2 —y 2 . Vj in Eq. ([T]) is the self-consistent one-electron potential, which 
is constructed using parameters of effective Coulomb interactions and the density matrix for 
the Mn 3d states. Generally, V« depends on both spin and orbital indices. 

In practice, the electronic low-energy model can be derived from the first-principles elec- 
tronic structure calculations, starting from the local-density approximation (LDA). 15 The 
construction of the model can be formulated rather rigorously in the basis of Wannier or- 
bitals for the Mn 3d bands. Then, is identified with the matrix elements of the LDA 
Hamiltonian in the Wannier basis. Thus, without Vj, the parameters Uj describe the LDA 
electronic structure for the Mn 3d bands. The parameters of effective Coulomb interac- 
tions for the Mn 3c? bands can be derived, also in the Wannier basis, using constrained 
random-phase approximation and/or the constrained LDA approach. For details, the reader 

6 



is referred to the review article (Ref. Il5l ). Then, the model can be solved in the mean-field 
HF approximation, which gives us the potentials Vj.— 

After the solution, the FE polarization can be obtained by applying the Berry-phase 
theory.— >ii Namely, the FE polarization is divided into the ionic (ion) and electronic (el) 
parts: 

p _ pion _|_ pel 

The ionic term reflects the non-cenrosymmetricity of the crystal structure itself and is as- 
sociated with the displacements (Ar^) of ionic charges (Zj) away from the centrosymmetric 
positions: 

pion = 7E^ Ar " ( 2 ) 

i 

where V is the primitive cell volume. The electronic term reflects the fact of the inver- 
sion symmetry breaking in the form of the wavefunctions, obtained from the solution of 
quantum-mechanical Schrodinger equations. It incorporates the effects of the magnetic in- 
version symmetry breaking and can take place even for centrosymmetric crystalline systems, 
provided that the inversion symmetry is broken by magnetic or some other electronic de- 
grees of freedom. The electronic term can be computed in the reciprocal space, by using the 
formula of King-Smith and Vanderbilt: 1 -- 

M 

J2 / (rik|V k |nk)dk, (3) 



>el ie 



n=l 



BZ 



where \nk) is the cell periodic wavefunction, the summation runs over the occupied bands 
(n), the k-space integration goes over the first Brillouin zone, and — e (e > 0) is the electron 
charge. In practical calculations, Eq. fl3]) is replaced by a discrete grid formula.— Eq. ([3]) 
can be also rewritten in terms of the Wannier function (w n ), constructed from |nk) in the 
real space:— 

M „ 

J r\w n (r)\dr. (4) 



M 

pel ' 



In all these equations, it is understood that P is the change of the polarization, obtained 
in the process of adiabatic lowering of the inversion symmetry— Moreover, the contribution 
of the low-energy bands (in our case, the Mn 3c? bands) is accounted by P . Therefore, the 
contribution of all other occupied states, which are not included to the low-energy model, 
should be described (at least, approximately) by P lon . Then, since the oxygen 2p band is 



fully occupied, it is reasonable to take Zq = — 2e, which corresponds to the formal valence 
state of 2 ~. On the other hand, all valence states of the rare-earth (RE) ions are empty. 
This should correspond to Zre — 3e. In the noncentrosymmetric P2inm structure, the 
Mn sites do not contribute to P lon .§ Therefore, the parameter Zy^ n is not important for our 
purposes. 

In the previous publications, this procedure was applied to the series of orthorhombic 
manganites. Particularly, the behavior of parameters of the low-energy model, derived from 



the first-principles electronic structure calculations, was discussed in Ref. 



116. 

such parameters for VMu() :f am be found in Supplemental Material of RefJs. The properties 
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of the magnetic ground state, obtained from the solution of the low-energy model in the H 



.0 



approximation, and corresponding behavior of the FE polarization were considered in Refs 
and 8. Note that a scaling factor was missing in the calculations of the FE polarization 
reported in Ref. U. This error was corrected in Ref. §]. 

As far as the FE polarization is concerned, the low-energy model reproduces results of 
the first-principles electronic structure calculations (Refs. |5|, |l7|, and ll8|) on a good semi- 
quantitative level. Moreover, the low-energy model was very helpful in clarifying details of 
the noncollinear magnetic ground state, which can be realized in orthorhombic manganites, 
namely: (i) the canting of spins and magnetic origin of the twofold periodic phase;- 1 ^ (ii) 
deformation of the spin-spiral texture, yielding FE activity in both two- and fourfold periodic 
systems;- (iii) the absence of the magnetic inversion symmetry breaking in systems with odd 
magnetic periodicity.- 

In this work, we will further rationalize the story by considering the DE limit for the FE 
polarization. 

Let us start with the ferromagnetic (FM) state, where each Vj is diagonal with respect 
to the spin indices, 



Vi 



vj 



x Vi- 
and Vj' 1 are the 5x5 matrices in the orbital subspace. The states with s =b are occupied 
by four electrons and the ones with s =1 are empty. Then, Vf can be identically presented 
in the form: Vf = A ox + AVf, where A ox is the intra-atomic exchange splitting between 
centers of gravity of the majority (f) and minority (\) spin states, and AVj describes the 
orbital splitting of unoccupied b-spin states. Moreover, four 3d electrons obey Hund's first 
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rule, which tend to form the state with the maximal spin S = 2. Therefore, besides on-site 
Coulomb repulsion (U), A ex will contain a large contribution, being proportional to the 
local magnetic moment (2S) and the intra-atomic exchange coupling (Jh)- This is the main 
reason why for many applications A ex can be treated as the largest physical parameter, and 
the DE limit corresponds to the extreme situation where A ex — > oo.— On the other hand, 
the splitting of unoccupied J,-spin states is considerably weaker. For example, in the HF 
approximation, it is caused by relatively small nonsphericity of the Coulomb potential. 

Therefore, when A cx — > oo, the details of (finite) splitting of the J,-spin states become 
unimportant and our first approximation is to replace AVf by Vj. It allows us to present 
Vi in the following form: 

H~v;+ f I , (5) 

\ A cx J 

where the orbital-dependent part (V.J) does not depend on the spin indices and the spin- 
dependent part does not depend on the orbital ones. Therefore, spin and orbital transfor- 
mations of Eq. flS]) can be treated separately. 

A typical example, illustrating the structure of the atomic 3d level splitting by the 
Coulomb and exchange potentials in the low-energy model, is shown in Fig. [TJ Typical 
values of A ex in manganites are about 4.5 eV, while the splitting of the 4--spin states is 
about 1.7 eV. The difference is not extremely large. However, as we will see in a moment, 
it is sufficient to justify the use of the DE limit for the FE polarization. 

As the next step, let us consider an arbitrary magnetic texture, where the directions 
of spin (e,) at each site of the lattice are specified by the combinations of polar (6^) and 
azimuthal (0,) angles: e; = (cos fa sin #j, sin fa sin cos 6^). Corresponding electronic struc- 
ture can be generated by the unitary transformation of Eq. using spin-rotation matrices: 

Vi -»> U(6i, fa)Vilfi(9i, fa), (6) 

where 

U(9i,fa) = 



cos f sin ^e~* 



— sin cos 2 



Here, it is assumed that the angles (di,(pi) are specified by magnetic interactions in the 
system (for the form of the optimized magnetic textures, the reader is referred to Refs. Q 
and |8|) and the one-electron potential for an arbitrary direction of spin can be obtained by 
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FIG. 1. Eigenvalues of the Hartree-Fock potential, as obtained in the 



ow-energy model for the fer- 
and |8| for the experimental 
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romagnetic phase of YMnC>3, HoMn03, and TbMnOs (results of Refs. 
Pbnm structure). A ex is the intra-atomic splitting between centers of gravity of the majority (t) 
and minority (|) spin states. 

using rigid spin rotations [Eq. fl6])] without additional self-consistency. This is a very good 
approximation in the case of manganites, because: 

(i) Due to the strong Hund's coupling, the local spin magnetization will always tend to 
stay in the saturated state. Therefore, the absolute value of this magnetization will 
only weakly depend on the direction of spins at other magnetic sites. 

(ii) The orbital configuration is rigidly fixed by the Jahn- Teller (JT) distortion and practi- 
cally does not depend on the type of the spin texture. For example, the energy splitting 
of the e g states, caused by the JT distortion, is about 1.5 eV, while typical strength of 
interatomic exchange interactions is of the order of several meV.— The exchange in- 
teractions can be additionally optimized by means of the orbital reconstruction, which 
works against the JT splitting.*^ However, the possible energy gain, caused by this 
reconstruction (typically, of the order of the exchange interactions themselves) is much 
smaller than the energy of the JT distortion. Thus, the orbital reconstruction does 
not occur. 
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The next step is to transform Eq. (J6]) to the local coordinate frame, corresponding to the z 
direction of magnetization at each site of the lattice. It leads to the following transformation 
of the transfer integrals: 

Then, taking the limit A cx — > oo, we obtain the well known DE model: 

which formulated in the subspace of the ^-spin states, in the local coordinate frame.— The 
prefactor £y is nothing but the 'f'N-ele 1116 ^ °f the product U^(8i,<f)i)U(8j,<f)j): 

£ a = cos — cos — + sin — sin ^i e ~^ l ^^\ 
s 3 2 2 2 2 

which satisfies the well known property: £y= 1 and for the ferromagnetically and antiferro- 
magnetically coupled spins, respectively. Therefore, in the DE limit, any antiferromagnetic 
(AFM) phase effectively breaks up into FM segments. For example, the description of the 
E-type AFM phase is reduced to the analysis of one-dimensional FM zigzag chains.- 1 ^ 

Next, we investigate abilities of the DE model for the description of the FE polarization. 
For these purposes, we calculate the electronic structure for the DE Hamiltonian [Eq. @], 
and then evaluate the electronic polarization, using the Berry-phase formula [the discrete 
analog of Eq. (E}].— ^ This procedure was applied to the series of orthorhombic manganites 
TbMnC>3, HoMnC>3, and YMnC>3 (and using both experimental and theoretically optimized 
crystal structure for the latter compound). 7 ' 8 The obtained polarization was compared with 
results of self-consistent HF calculations for the same low-energy model, but without ad- 
ditional approximations associated with the use of the DE limit. Typical results of such 
calculations are illustrated in Fig. |2] for the Pbnm phase of YMnC>3 (other systems show 
very similar behavior). More specifically, we consider a twofold periodic magnetic texture, 
which is explained in Fig. EJ^c), and keep the AFM coupling between adjacent planes z = 
and z = c/2, as explained in Fig. Ufa). Then, = and 180° correspond to the AFM 
alignment of the E-type, while = 90° corresponds to the spin-spiral alignment. For this 
geometry, the FE polarization should be parallel to the orthorhombic a axis.- In the DE 
model itself, we consider two levels of approximations. In the first case (denoted as 'DE 
LDA'), we neglect Vj and consider only the crystal-field splitting and transfer integrals, 
derived from the LDA band structure. Then, the transfer integrals are modulated by 
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FIG. 2. (Color online) (a) E-type antiferromagnetic texture, (b) Behavior of electronic polarization 
in YM11O3 upon rotation of magnetic moments as obtained in the self-consistent mean-field Hartree- 
Fock approximation (total HF); in the double exchange model for the LDA band structure (DE 
LDA); and in the double exchange model with the Hartree-Fock potential vj (DE +17). In the 
rotated texture, the directions of spins at the sites 1 and 3 were fixed, while the spins at the sites 
2 and 4 were rotated by the angle <j), as explained in panel (c). The planes z = and z = c/2 were 
coupled antiferromagnetically. 

as requested by the DE model. In the second case, we consider the full DE Hamiltonian, 
Eq. ©, including Vj (denoted as 'DE +[/'). All magnetic solutions are insulating. There- 
fore, we can use the Berry-phase formula for the analysis of P cl . The DE LDA scheme 
overestimates the electronic polarization by about 50 %. Nevertheless, this is to be ex- 
pected, because LDA underestimates the band gap. Therefore, the FE polarization should 
be generally larger. Similar behavior was found in the first-principles calculations.- 1 ^ The 
analytical expression, explaining the band-gap dependence of P el , will be derived in the next 
section. The band-gap problem is corrected by Vj. Therefore, the FE polarization, derived 
in the DE +U scheme, is smaller. Moreover, results of self-consistent HF calculations for 
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the electronic polarization are well reproduced by the DE +U scheme: although P in the 
approximate DE +U scheme is systematically smaller, the typical difference, which was 
obtained for all considered systems, is less than 15 %. 

This is our main observation and also the main motivation of the rest of our work. By 
considering the DE limit, we will slightly lose in the accuracy. But instead we will be able 
to rationalize the problem and derive several analytical expressions for the FE polarization 
in orthorhombic manganites. Our analysis will also clarify results of the low-energy model 
and first-principles calculations. 

III. RESULTS 

We start with the analysis of the E-type AFM phase. As was pointed out above, in 
the DE limit, the FE AFM E-pha.se breaks up into one-dimensional FM zigzag chains. 
Therefore, the key moment for understanding the origin of the FE activity in the E'-phase 
is the analysis of isolated zigzag chain.-^ In Sec. IIII A\ we start such an analysis with the 
simplest but analytically solvable model for the e g electrons. In Sec. IIII Bl we will derive 
parameters of such a model, starting from a more general five-orbital model, which was 
obtained from the first-principles calculations.- 1 ^^ From the analysis of this model we will 
conclude that the situation, realized in most of the electronic structure calculations (even in 
ordinary LDA), corresponds to the limit of large energy splitting A between atomic e g states, 
which incorporates the effects of the JT distortion and (optionally) the on-site Coulomb 
repulsion. Then, by considering the large-A limit, in Sec. IIII CI we will derive an analytical 
expression for the FE polarization, which is based on the five-orbital model. This expression 
explains the functional dependence of P el on the relative directions of spins and the form 
of nearest-neighbor transfer integrals. In Sec. IIII Dl we will analyze relative directions of 
electronic and ionic polarizations in the noncentrosymmetric P2inm structure and point 
out on the problem of structural optimization, which apparently exists in some of the first- 
principles calculations, where the directions of noncentrosymmetric atomic displacements 
are inconsistent with the type of the orbital ordering, realized in the FM zigzag chain. 
In Sec. IIII El we discuss the possibility of switching the FE polarization by changing the 
magnetic texture: we argue that, even in the twofold periodic texture, there is another type 
of the AFM zigzag ordering, which leads to a finite FE polarization along the orthorhombic 
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FIG. 3. (Color online) Geometry of the zigzag chain for the square lattice and the occupied e g 
orbitals of the 3x 2 —r 2 and 3y 2 —r 2 type. Cubic and orthorhombic axes are denoted as xy and ab, 
respectively. 

c axis. However, the value of this polarization is expected to be small. 

A. Analytically solvable model for the e g electrons in the zigzag chain 

The zigzag chain consists of the two groups of sites: the lower corner sites 1 and the 
upper corner sites 2 (see Fig. [3]). The orthorhombic translation a transforms each group 
to itself (the translated sites are denoted as 1' and 2', respectively). It is assumed that 
the lattice distortion stabilizes some e g orbitals at the sites 1 and 2, which will be denoted 
as |l)i and 1 1) 2 , respectively. The orthogonal to them e g orbitals are denoted as \2)i and 
|2)2, respectively. Furthermore, it is assumed that there is a symmetry operation (S), which 
transforms the zigzag chain to itself and which consists of the 180° rotation around the a 
axis (C 2 ) with consequent translation. S will transform site 1 to site 2, and vice versa. 
For the Phnm structure (and with some appropriate choice of the origin), such symmetry 
operation is {C 2 \a/2+b/2} (where the first part stands for the rotation, and the second part 
specifies the translation) while for the P2inm structure, it is {C 2 \a/2+c/2}. It is important 
that both symmetry operations include the translation a/2. Then, it is convenient to work 
in the local basis, corresponding to the diagonal presentation of the e g level splitting, such 
that S would transform the basis functions of the site 1 to the ones of the site 2, and vice 
versa. Our idea is that, although we have two different sites, with such choice of the basis 
functions, the Hamiltonian becomes periodic with the period a/2 and the problem can be 
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treated as if it would have only one site in the primitive cell. Similar idea was used for the 
analysis of the CE AFM state in the half-doped manganites.— Rather generally, these basis 
functions can be chosen in the form: 

= -cos/3|32 2 -r 2 ) 1 -sm/3\x 2 - y 2 )i, (8) 

|2)i= sin/3|32 2 -r 2 ) l - cos f3\x 2 -y 2 ) 1 (9) 

at the site 1, and 

|1) 2 = -cos/3|3^ 2 -r 2 ) 2 + sin/3|x 2 -?/ 2 ) 2 , (10) 

|2) 2 = sin/3|3z 2 -r 2 ) 2 + cos/3|x 2 -i/ 2 ) 2 (11) 

at the site 2, where — 7r/2 < f3 < ir/2. \(3\ = 60° corresponds to the ideal square lattice, 
subjected to the JT distortion. Here, it is assumed that the direction of this distortion 
is determined by anharmonic electron-lattice interactions, which stabilize orbitals of the 
type |1) with \/3\ close to 60°.— >2i Then, all deformations of the orbital ordering pattern are 
described by the single parameter /3. Here, we continue to use the notations |3z 2 — r 2 ) and 
\x 2 —y 2 ) for the e g orbitals, although it should be understood that they are valid only for 
the ideal square lattice, and more generally we have in mind some |3z 2 — r 2 )-like orbitals, 
which transform to each other as S\3z 2 — r 2 )i = |3^ 2 — r 2 ) 2 , and some \x 2 — y 2 )-like orbitals, 
which transform to each other as S\x 2 — y 2 )\ = — \x 2 — y 2 )i- Such deformations of the ideal 
e g orbitals can be caused, for example, by buckling distortions. It is easy to check that 
(3 = -60° yields |1>! = |3x 2 -r 2 )!, \2) x = ^-z 2 ),, |1) 2 = \3y 2 -r 2 ) 2 , and |2) 2 = \x 2 -z 2 ) i; 
while p = 60° yields |l)i = |3y 2 -r 2 )i, |2)j = |^ 2 -x 2 )i, |1) 2 = |3x 2 -r 2 ) 2 , and |2) 2 = 
\z 2 —y 2 ) 2 - Then, although in realistic situations, |/3| can deviate from 60°, we will say that 
/3<0 corresponds to the 3x 2 — r 2 /3y 2 — r 2 type of the orbital ordering (referring to the type of 
the occupied orbitals at the sites 1/2), while /3>0 corresponds to the 3y 2 —r 2 /3x 2 —r 2 type 
of the orbital ordering. 

As for the transfer integrals between e g orbitals, we again consider a more general case 
and write them in the following form: 



for the bond 1-2', and 



'12' 



12 



-- (l- | sinp\a x - cos/35-^ , (12) 
-^(i+|sin/3|<7 x -cos/?<7*), (13) 
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for the bond 1-2, in terms of the pseudospin Pauli matrices a x , a y , and a z , and the 2x2 
identity matrix I. Throughout this section, all energies are in the units of two-center integral 
to of the dda type.— The form of t^ is suggested by the dda transfer integrals in the ideal 
square lattice, which again corresponds to = 60°. Therefore, it is assumed that all 
deviations from the ideal square lattice are described by the single parameter (3, similar to 
the orbital ordering. Note also that Eqs. ( IT2|) and (1131) satisfy the idempotency condition 
(tij) 2 = Uj, which holds for the dda transfer integrals in the square lattice. 

Thus, in our model, the orbital ordering and the transfer integrals are described by 
the same parameter /3. Generally speaking, these are different quantities, which should be 
specified by two different sets of parameters. Nevertheless, in the analytical model, one would 
always like to reduce the number of independent parameters to the minimum. Moreover, the 
use of the single parameter /3 is indeed a very reasonable approximation for our purposes: 

(i) At least for the ideal square lattice, the orbital ordering and the transfer integrals can 
be described by the same \/3\ = 60°. Thus, there is the reference point where our 
construction is exact; 

(ii) Small deviations from the ideal case are treated as an approximation and we have 
some freedom to decide the form of this approximation. In Sec. IIIIBI we will show 
that typical deviations of from 60° are not large and, therefore, our approximation 
is robust; 

(iii) According to Eqs. (j!2p and (fl~3j) . the transfer integrals do not depend on the sign of /3 
(although the orbital ordering does). This is the very important requirement, because 
the phases of transfer integrals is determined solely by the geometry of the zigzag 
chain and should not depend on the type of the orbital ordering. 

After the transformation to the local basis, given by Eqs. (l8])- ffTTj) . the transfer integrals 
become: 

ti2' = t 2 i = - ^cos/31 + sin 2/3(72; — i\ sin f3\a y — cos2/3a z ^j , (14) 

where t-2'i = ti2 = t^. Thus, the transfer integrals are indeed periodic with the period a/2 
and, in the reciprocal space, the problem is reduced to the analysis of the 2x2 Hamiltonian 
of the form: 

H{k) = e(Jfe) + d(Jfe) ■ <r, 
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where e(k) = cos/3 cos(£;a/2), and components of the vector d = (d x ,d y ,d z ) are given by 
d x = sin 2/3 cos(fca/2), d y = | sin/3| sin(fca/2), and d z = — cos 2/3 cos(/ca/2) — A/2. The 
parameter A in d z is the intra-atomic energy splitting between e g states, caused by lattice 
distortions and Coulomb interactions. This result can be also viewed as if the transformation 
(18j)- (lTTj) would "straightened" the zigzag chain and made it equivalent to a linear chain, 
but with different transfer integrals operating in the positive and negative directions of a. 
Because of the condition ti2 = t^, the transfer integrals are generally not centrosymmetric 
with respect to the atomic sites and the system will develop a finite electronic polarization. 
Nevertheless, in the limit A— >oo, the basis orbitals of the type '2' are projected out. Then, 
the transfer integrals between orbitals of the same type '1' are just scalars, and the condition 
ti2 = % becomes equivalent to ti 2 = t 2 i- Thus, in the limit A— >oo, the problem should 
become centrosymmetric. From this point of view, it is logical to consider the limit A— >oo 
as the reference point for the electronic polarization. 

The eigenvalues of 'H(k) are given by E±(k) = e(k)±\d(k)\, and the eigenvector, corre- 
sponding to the lowest occupied band, satisfies the condition: [d(/c) ■ <y+|d(A;)|]| — , k) = 0.— 
Then, |— , k) can be taken in the form: 



and 7(/c) = arctan^/da,). 

At the half-filling (one e g electron per each Mn site), the zigzag chain is a band insulator. 
This property holds even for A = due to specific form of the dda transfer integrals. 24 
Moreover, the reciprocal lattice vector of the "straightened" chain is G = 4n/a, and |— ,k) 
is a periodic function of G. Therefore, the electronic polarization can be computed directly, 
using the formula of King-Smith and Vanderbilt ~i Note that in this section, it is more con- 
venient to work with the electric dipole moment, rather than with the polarization density. 
Therefore, Eq. fl3]) was additionally multiplied by the primitive cell volume V. Neverthe- 
less, unless it is specified otherwise, we will use the same notations for this quantity and 




where 
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continue to call it "the polarization". Then, we obtain the following expression for the FE 
polarization parallel to the orthorhombic a axis (per two Mn sites in the zigzag chain): 



E 



which can be further transformed to 



4Wo |d(*)|[|d(*)| -<(,(*)] 

where the subscript E means that this polarization corresponds to the E-type AFM phase 
in the DE limit. 

Thus, we immediately recognize that when the orbital ordering changes from 3x 2 —r 2 /3y 2 —r 2 
(/3<0) to 3y 2 —r 2 /3x 2 —r 2 ((3>0), the polarization changes its sign. 
Then, it is straightforward to find that 

lim Pi = ^™ 
a^o+ E sin/3 2 



and, therefore, |P|J| = ea/2 (see Ref. l25l ). Then, since P§ is well defined modulo ea,— the 
values of P§ and — P§ for A = are equivalent. Such a situation means that the system 
possesses the inversion symmetry, but the inversion centers are located in the middles of the 
bonds— Thus, by removing the JT distortion from our model, we effectively create a new 
inversion center. This is indeed the case for the model considered above: since tu = £21 and 
tw = tyi [see Eqs. (Tl2"]) - (113p ]. the transfer integrals are centrosymmetric with respects to 
the middles of the bonds. 

In the limit A — > 00, we have 

^ ea|sM|sin2fl 
A 2 

This result also has a transparent physical meaning and can be easily understood by starting 
from the expression 

pf = -2e / xw 2 {x)dx, (17) 



in terms of the Wannier functions,-^ where the prefactor '2' stands for the number of Mn 
sites in the primitive cell of the zigzag chain. Let us consider the limit A — > 00, where 
l^oo) = and it is centered at the site 1 (see Fig. |3]). Then, in the first order of 1/A, this 
Wannier function will have a finite tail, spreading to the neighboring sites 2 and 2', which 
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are located at x = —a/2 and a/2, respectively. In the first order of perturbation theory, this 
tail is proportional to the transfer integrals [Eq. ffT4j) ] from the occupied orbital |l)i to the 
subspace of unoccupied orbitals |2) at the sites 2 and 2'. Then, by assuming that all weights 
of w 2 (x) are accumulated at the lattice points (that is the meaning of the "lattice model"), 
one can write that 



are the weights of w 2 (x) at the sites 2 and 2'. By substituting this w 2 (x) into Eq. ffTTj) . 
we again arrive at Eq. ( fT6l) . Thus, in terms of these arguments, the polarization is finite 
because q + ^ g_. Alternatively, one can say that due to the asymmetric electron transfer, 
the Wannier centers are shifted from the centrosymmetric atomic positions.™ For a given 
A, the difference (q + — qJ) depends on the value of and takes the maximal value when 
= |/3 m ax| = arctan \/2 (about 54.7°). 

The behavior of electronic polarization as the function of intra-atomic energy splitting 
between e g states is summarized in Fig. HI 

A very similar model of the FE polarization in orthorhombic manganites was considered 
by Barone et al.— The advantage of our approach is that we were able to reduce the problem 
to the 2x2 Hamiltonian in the reciprocal space and to solve it analytically. Such an analysis 
provides a transparent physical picture for the behavior of the FE polarization. Therefore, 
we would like to stress briefly the difference between our results and the ones by Barone et 
al. First, the behavior of polarization, obtained by Barone et al., is very different from ours: 
it is zero for A = and approaches ±ea/2 for A — > oo. Nevertheless, such a difference can 
be easily understood by the different choice of the reference point in the calculations of Pjh 
A — > oo in our work and A = in the work of Barone et al. Another discrepancy is related 
to the functional dependence of the orbital ordering and the FE polarization on A: in the 
work of Barone et al, these two quantities become finite starting only from some critical 
value of A. We believe that such a behavior is counterintuitive (at least, in the framework 
of the considered model) and the orbital ordering, as well as the FE polarization, should 
evolve continuously starting from A = (see also the analysis of the orbital ordering for 
similar model, reported in Ref. uM . 




w 2 {x) = (1 - g_ - q+)8{x) + q-5(x + a/2) + q + 5{x - a/2), 



where 
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A/|sin2/? sin/? | (units of t Q ) 

FIG. 4. (Color online) Electronic polarization (more precisely - the electric dipole moment) for the 
isolated zigzag chain as the function of intra-atomic energy splitting between e g states. Upper part 
corresponds to the 3y 2 — r 2 /3x 2 — r 2 type of the orbital ordering (/3>0) and lower part corresponds 
to the 3x 2 — r 2 /3y 2 — r 2 type of the orbital ordering (/3<0). The values obtained for YMn03 are 
shown by open symbols. The points A, C, E, G, I denote the bare LDA values, obtained for the 
experimental P2\nm and Pbnm structures, and three theoretical structures, obtained in LSDA 
and LDA+J7 with U = 2.2 and 6.0 eV, respectively. Similar values, obtained after adding the HF 
potentials, are denoted as B, D, G, H, and J, respectively. The values obtained for the experimental 
Pbnm structure of HoMnOs are shown by filled symbols: the point K denotes the bare LDA value 
and the point L takes into account the effect of the HF potential. 

B. Parameters of the e g model and values of electronic polarization for YMnC»3 
and HoMn0 3 

In this section, we evaluate parameters of the e g model for realistic compounds, such as 
YMnC>3 and HoMnOs. For these purposes, we do the following: 

(i) We start with the realistic low-energy model, derived for the Mn 3d bands of YMnOs 
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and HoMnOs on the basis of first-principles electronic structure calculations (results 



of Refs. 



16 and 



(ii) Then, we pick up parameters of the model for the single zigzag chain, propagating along 
the orthorhombic a axis (and assuming that, in the DE limit, all transfer integrals in 
the directions b and c are blocked by the E-type AFM ordering); 

(iii) Solve the electronic structure problem for the isolated zigzag chain; find eigenvalues 
and eigenfunctions; 

(iv) Construct the Wannier functions for the upper lying e g bands. For these purposes, we 
use the projector-operator technique and trial orbitals, obtained from the diagonaliza- 
tion of the site-diagonal part of the density matrix;— 

(v) Find parameters of the e g model in the obtained Wannier basis; 

(vi) Transform the parameters to the crystal-field representation, which diagonalizes the 
site-diagonal part of the e g model; 

(vii) Fit the transfer integrals for the bond 1-2' in terms of t and (3, by using the functional 
dependence given by Eq. flSJ). Meanwhile, the splitting A between the e g states is 
obtained from the site-diagonal part. 

For YMnC>3, we have considered several crystal structures, which were previously dis- 
cussed in Ref. 



(i) The experimental Pbnm and P2\nm structures, reported in Ref. [17 



(ii) Three theoretical P2inm structures, which were optimized in the local-spin-density 
approximation (LSDA) and LDA+U with U = 2.2 and 6.0 eV by assuming the collinear 
E-type AFM alignment without SO interaction. The results of this optimization can 
be found in Ref. 



For HoMnOs, we use the experimental Pbnm structure, reported in Ref. 127 . 

Parameters of the e g model, obtained from the fitting, are summarized in Table HI Eq. (I14p 
captures main details of transfer integrals between the nearest neighbors. The largest devi- 
ation from the ideal = 60° case was found if one uses the bare LDA parameters, derived 
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TABLE I. Parameters of the e g model for the isolated zigzag chain, derived for H0M11O3 (HMO) and 
different structures of YMnC>3 (YMO): the experimental Pbnm and P2\nm structures, reported 



in Ref. 1 171. and three P2\nm structures, which were theoretically optimized in LSDA and LDA+£7 
with U= 2.2 and 6.0 eV by assuming the collinear E'-type AFM alignment (results of Ref. 8). In 
this Table, io is the effective two-center integral ddcr, A is the intra-atomic splitting between e g 
states, and /3 specifies the form of the transfer integrals in the Mn-Mn bonds. The values, obtained 
by using bare LDA parameters are denoted as 'LDA', and the ones after adding the Hartree-Fock 
potential are denoted as l +U\ 





t (meV) 


A 


(eV) 


P (de* 


;rees) 


LDA 


+u 


LDA 


+u 


LDA 


+u 


HMO (Pbnm, Exp.) 


341 


353 


1.52 


2.15 


-54.0 


-55.3 


YMO (Pbnm, Exp.) 


335 


348 


1.53 


2.18 


-54.2 


-55.8 


YMO (P2xnm, Exp.) 


334 


346 


1.54 


2.15 


-54.1 


-55.7 


YMO (P2 1 nm, LSDA) 


405 


412 


0.92 


1.95 


57.6 


59.1 


YMO (P2mro, U= 2.2 eV) 


361 


370 


1.37 


2.41 


55.1 


57.3 


YMO (P2mm, U= 6.0 eV) 


348 


359 


1.30 


2.04 


54.4 


56.1 



for the experimental Pbnm structure of HoMnCV In this case, the agreement between the 
original matrices i±2> and results of the fitting using Eq. (JI4"]) is the worst: 

/ 119 -328 \ / 153 -300 

ti2' = and 

V —12 39 I \ -24 47 

before and after the fitting, respectively, in units of meV. On the other hand, /3 becomes close 
to 60° if one uses theoretical LSDA crystal structure of YMnC>3 and takes into account the 
additional level splitting, caused by the HF potential. In this case, the agreement between 
the original and fitted matrices is nearly perfect. Nevertheless, we would like to emphasize 
that the analytical expression, given by Eq. (ITo"]) . with the parameters, derived from the 
fitting, excellently reproduces the behavior of electronic polarization, obtained in the same 
e g model but without fitting (see Fig. [5]). Thus, deviations of transfer integrals from Eq. (fl4"]) 
are relatively unimportant for the analysis of the FE polarization. 
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FIG. 5. Electronic polarization (more precisely - the electric dipole moment) in the e g model for 
the isolated zigzag chain as obtained by using bare LDA parameters for H0M11O3 (solid line) and 
after the parametrization of transfer integrals using Eq. (|14p . The vertical line shows the bare LDA 
value of A/to- The parameters are taken from Tabled 

It is interesting to note that \P§\ = ea/2 when A — > + (see Fig. [5]). This means that in 
the limit A — > 0, the system behaves such as if it would be centrosymmetric with respect 
to the bond centers,— even despite the fact that the space groups Pbnm and P2inm (so as 
the transfer integrals) do not have such symmetry. Apparently, such a behavior is related 
to a more general symmetry of the transfer integrals. 

Then, we take the values of FE polarization, obtained for the e g band (without fitting) and 
also plot them on Fig. HJ As for the abscissa coordinates, we use results of Table HI We can 
clearly see that all these values fall on the analytical dependence, derived for the e g model. 
The main parameter, which controls the value of the FE polarization, is the ratio A/t . 
The /3-dependence of P§ is less important. This result is very natural and will be discussed 
in a moment. Moreover, the physically relevant situation, realized in the orthorhombic 
manganites, always corresponds to the limit of large A. This is another important finding, 
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which will allow us to further rationalize the behavior of the FE polarization in Sec. II II CI 
The polarization has different sign for the experimental and theoretical structures, that 
indicates at different types of the orbital ordering in the zigzag chain. In the Pbnm phase, 
all zigzag chains are equivalent, and in Fig. H]we simply picked up the one with the same 
orbital ordering as in the P2inm phase. However, in the P2inm phase, the type of the 
zigzag chain is uniquely defined (as the one with larger Mn-Mn distances, which stabilize 
the FM coupling in the zigzag chain). Therefore, the sign difference between experimental 
and theoretical values of PfJ in the P2 1 ram phase indicates at a serious problem, which 
may exist in the first-principles calculations. The problem will be discussed in details in 
Sec. lfflDl 

Then, all values of \/3\ are close to |/3 max | ~ 54.7°, which corresponds to the maximum 
of \P§\ (see Tabled]). Therefore, any deviation of P§((3) from P§((3 max ) will be only of the 
order of ((3— /3 max ) 2 - Thus, all the effects of (3 on PjJ will be small. This can be clearly 
seen in Fig. El where we plot P§ versus A/to, using different sets of parameters for the e g 
model and varying A: all lines, corresponding to different crystal structures and different 
levels of approximation for the on-site interactions (with and without the HF potential), are 
practically undistinguishable. This means that, in reality, P§ is controlled by only two sets 
of parameters: (i) the ratio A/to, and (ii) the lattice parameters a, b, and c, which determine 
the value of the scaling factor a/V in the polarization density. The /3-dependence of P§ is 
relatively unimportant. 

From the physical point of view, the /3-dependence of the transfer integrals is related to 
the buckling of the Mn-O-Mn bonds. Then, the above result suggests that P§ does not 
explicitly depend on the Mn-O-Mn angles: the latter can contribute to PjJ, but only via 
other model parameters (such as to), which depend on these angles. This finding is consistent 



with the conclusion of Ref. 
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based on the first-principles electronic structure calculations. 
Finally, we briefly explain the correspondence between the values of the electric dipole 
moment in Fig. H] and the polarization density. Let us consider the experimental Pbnm 
structure of YMnC>3. Then, the value — 0.022ea, which takes into account the effect of the 
HF potential, corresponds to the polarization density of about —1.65 /zC/cm 2 . It should be 
remembered that it is only the contribution of the e g band alone. In order to obtain the 
total polarization for the five-orbital model, it should be combined with the contribution 
of the t 2g band. This yields the total polarization —0.84 /zC/cm 2 , which agrees with the 
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A (units of ? ) 



FIG. 6. (Color online) Electronic polarization (more precisely - the electric dipole moment) versus 
A/to, as obtained using various sets of parameters for the e g model. Results of bare LDA and 
after including the Hartree-Fock potential are shown by filled and open symbols, respectively. The 
positions of A/to f° r different systems are shown by capital letters. The points A, C, E, G, I 
stand for the bare LDA values of A/to, corresponding to the experimental P2\nm and Pbnm 
structures, and three theoretical structures, obtained in LSDA and LDA+C7 with U = 2.2 and 6.0 
eV, respectively. Similar points, obtained after adding the Hartree-Fock potential, are denoted as 
B, D, G, H, and J, respectively. The points K and L correspond to the Pbnm structure of HoMn03, 
obtained in the bare LDA and after including the HF potential, respectively. 



value for the E'-type AFM state (for <fi = 180°) in Fig. [2J Thus, the contributions of the 
tig and e g bands have opposite sign and partially cancel each other, in agreement with the 
first-principles calculations.— In the rest of this work, we will deal with the total polarization 
density, including the effect of both £ 2 g and e g bands. 
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FIG. 7. (Color online) General twofold periodic magnetic texture in the ab plane of orthorhombic 
manganites, which remains invariant under the symmetry operation S = {C^\a/2+b/2} of the 
space group Pbnm. Solid and dotted lines denote two types of magnetically inequivalent bonds. 

C. Electronic polarization in the five-orbital model: simple analytical expression 

Now, we will generalize results of two previous sections and derive an approximate, but 
very transparent expression for the electronic polarization in orthorhombic manganites with 
a general twofold periodic magnetic texture. Our starting point is that the behavior of 
electronic polarization in realistic compounds corresponds to the limit of large A. This limit 
can be justified even without on-site Coulomb interactions (i.e., considering the ratio of 
transfer integrals to the crystal- field splitting in bare LDA), and is additionally strengthened 
after including the Coulomb interactions. Thus, the central quantity, which we should 
evaluate in the second order of 1/A, is the weight wf^j, transferred from the Wannier 
orbital at the site i to the neighboring site j. Moreover, since electronic polarization is equal 
to zero for the fully occupied band, it is more convenient to start with the unoccupied e g 
orbitals and consider the transfer integrals to the subspace of three t 2g and one e g occupied 
orbitals at each of the neighboring sites. This procedure should give us — P el . 

The transfer integrals obey certain symmetry rules and, in the DE model, are additionally 
modulated by £y. More specifically, we consider a planar magnetic texture which is shown in 
Fig. [7J The periodicity of this texture along the orthorhombic axes is a and 2b, respectively. 
The directions of spins are specified by three azimuthal angles: 02, 03, and 04 (while X = 
is treated as the reference point). Moreover, we assume that the DE Hamiltonian remains 
invariant under the symmetry operation S = {C^\a / 2+b / 2} , which transforms the bond 1-2 

26 



to 4'-3, the bond 3-2 to 4-5, etc. In the DE model, it imposes additional conditions on the 
azimuthal angles: cos 4f = cos ^ 3 ~^ 4 and cos 4f = cos ^ 3 ~^ 2 , which are satisfied if 3 = 02±04 
(modulo 27r). Thus, the magnetic texture is specified by only two independent parameters 
02 and 04, similar to the magnetic texture obtained in the mean-field HF calculations with 
the SO coupling.-^ 

Then, we consider some central site (say, site 3 in Fig. [7J and evaluate its contribution to 
the electronic polarization, which is caused by the Wannier weight transfer to the neighboring 
sites 4', 4, 2, and 2', located at (a+b)/2, — (a— b)/2, — (a+6)/2, and — (a— 6)/2, respectively. 
In the second order of 1/A (and apart from the proportionality coefficient, which will be 
specified later), the contribution of the site 3 to the vector of electronic polarization can be 
written as 

P? ~ 6 - cos 2 ^ [(a + b)w%_> 4 , -(a- b)w^ 4 ] + 

\ cos 2 y [(a - b)wj^ 2 , - (a + b)wj^ 2 ] , (18) 

where wf^j is proportional to the sum of squares of the transfer integrals from the 
unoccupied orbital 5 at the site i to the occupied orbitals 1-4 at the site j: wf^j = 
[(ty) 2 + (tf/) 2 + (tff) 2 + (tf/) 2 ] /A 2 . These transfer integrals should be calculated in the 
'crystal-field representation', that diagonalizes the site-diagonal part of the one-electron 
Hamiltonian. The parameter A is understood as the energy difference between the unoccu- 
pied orbital 5 and the center of gravity of occupied orbitals 1-4 (see Fig. [8]). Thus, in this 
analysis, we neglect the splitting between the occupied orbitals, which is smaller than A. 
Then, in the Pbnm structure, each Mn site is located in the inversion center. Therefore, 
wf_>,- in the bonds 3-4' and 3-2 (as well as 3-2' and 3-4) are equivalent, and Eq. ( {TBI can be 
further transformed to 

P^ 1 ~ | (cos 02 - cos 4 ) [(a + b)wl^ A , - (a - b)w 2 3 _, 4 ] . (19) 

Similar analysis can be performed for another Mn site in the primitive cell (say, site 4' 
in Fig. [7]). Moreover, since the sites 3 and 4' are connected by the symmetry operation 
S = {C 2 |a/2+b/2}, using Eq. ffl9l) . one can immediately obtain that 

P4J ~ ^ (cos 02 - cos 4 ) [(a - 6)^3^4/ - (a + &)u> 2 ^ 4 ] . 

Then, the total polarization P cl = 2(P3 1 +P2!) can be evaluated as 

P el = y (cos 02 - cos0 4 ) [wl__> 4 , - wl_> 4 ] a. 
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FIG. 8. Splitting of 3d levels for the experimental Pbnm phase of YMnO"3. The values, obtained 
using bare LDA parameters of the low-energy model are denoted as 'LDA', and the ones obtained 
after adding the Hartree-Fock potential are denoted as '+£/'. A is the energy splitting between 
the unoccupied orbital 5 and the center of gravity of occupied orbitals 1-4. 

Here, V is the primitive cell volume, containing four Mn sites (two in each of the ab planes, 
located at z— and c/2, which is reflected in the additional prefactor 2 in the above expres- 
sion). Finally, by applying the symmetry operation S = {C%\a/2+b/2}, the sites 3 and 4 
can be transformed to the sites 4' and 3, respectively. Thus, P el can be expressed through 
the transfer integrals in only one nearest-neighbor (NN) bond 3-4' (or in any equivalent to 
it bond): 

P cl = ^(cos0 2 -cos0 4 )Pl (20) 

where 

pi = | - « (2i) 

is the electronic polarization in the £7-type AFM state. For an arbitrary direction of spin at 
the site 1, the angular dependence (cos0 2 — cos0 4 ) in Eq. f[2"Ul) should be replaced by a more 
general expressions ei • (e 2 — e 4 ). Eqs. ( 1201) and ( 1211 allow us to rationalize many aspects of 
the multiferroic activity in manganites with the twofold periodic magnetic texture, namely: 

(i) P el is parallel to the orthorhombic a axis; 
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(ii) If 04 = 02+tt, P el is proportional to cos 02, which nicely explains the functional de- 
pendence of P el (0) in Fig. Et^b) and in the first-principles calculations for the same 
magnetic geometry (Ref. |5|); 

(iii) P el vanishes in the homogeneous spin-spiral state (02 = 7r/2 and 04 = 37r/2). This is a 
very natural result from the viewpoint of the DE physics: in the spin-spiral texture, all 

are the same. Therefore, all bonds remain equivalent, and the inversion symmetry 
of the DE Hamiltonian is not broken; 

(iv) Since t,j = t£, P*| can be also presented in the form 



PJ (22) 



where (v + , vJ) is the scalar product of the 4-dimensional vectors v± = (v±, v\, v\, v±), 
constructed from symmetric (+) and antisymmetric (— ) parts of the transfer integrals: 
v™ = tf™±t^ 5 . Thus, in order to have finite P el , the matrix of transfer integrals should 
have both symmetric and antisymmetric components. 

Let us evaluate Pg = (P|,0,0), using Eq. ( 1221) . for the experimental Pbnm phase of 
YMnOs. In this case, the unit cell volume is V = 224.13 A 3 and the orthorhombic lattice 
parameter is a = 5.245 A.— Then, for the bare LDA band structure, we have: A = 2.40 
eV (see Fig. EJ, v + = (-125,18,15,336) meV, and v_ = (99,-49,-25,-314) meV (all 
parameters of the low-energy model for YMnC>3 can be found in Supplemental Material 
of Ref. l8|). By substituting all these values in Eq. ( 122]) . we obtain P§ = —1.55 /iC/cm 2 , 
which agrees very well with the value of —1.53 /xC/cm 2 , obtained directly from the Berry- 
phase formula [Eq. (jHJ)], without additional approximations (apart from the DE limit). For 
the more realistic case, including the effect of the HF potential, we have: A = 3.45 eV, 
v + = (6,-117,26,335) meV, and v. = (10,91,-24,-319) meV. Then, Eq. ([22]) yields 
PjJ = —0.74 /iC/cm 2 , which is again consistent with the value of —0.85 /iC/cm 2 , obtained 
directly from the Berry phase formula [Eq. Q]. Moreover, the values of the scalar product 
(v + , V-) appear to be very close when they are calculated with and without the HF potential: 
—0.118 and —0.119 eV 2 , respectively. This result is very natural because the form of the 
crystal-field orbitals in orthorhombic manganites is mainly controlled by the JT distortion: 
the latter is large and thus 'decides' which orbitals will be occupied and which will not. On 
the other hand, the effect of on-site Coulomb interactions, being inversely proportional to 
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U,— is considerably smaller. Thus, although the Coulomb interactions contribute to the 
splitting between occupied and empty states (see Fig. [8]), they practically do not change 
the subspace of occupied orbitals. Therefore, the construction (v + ,v_), which is evaluated 
in the crystal-field representation, will not strongly depend on whether it is calculated with 
or without the HF potential. In such a situation, the absolute value of P§ will be mainly 
controlled by the parameter A in the denominator of Eq. (|22|) . 

Furthermore, A can be presented in the form: A = Ajt + Au, where Ajt and Au take 
into account the effects of the bare JT distortion and the on-site Coulomb interactions, 
respectively. In the example considered above, Ajt is the LDA level splitting and Au is 
the additional splitting, caused by the HF potential (see Fig. [8]). Then, if P|(0) is the FE 
polarization in LDA, the effect of on-site Coulomb interactions on P|J can be evaluated using 
the following scaling relation: 

P|(A [/ ) = P|(0)/(1 + A l/ /A JT ) 2 , 

which was observed in many LDA+[7 calculations, treating the on-site Coulomb repulsion 
U as an adjustable parameter.-*^ 

Finally, it is instructive to evaluate P cl = (P el , 0, 0) for the noncollinear magnetic ground 
state of YMn0 3 using Eq. ( I2"0j) . This magnetic ground state was obtained in Ref. |§] by 
solving mean-field HF equations with the relativistic SO interaction. For the Pbnm phase 
of YMn0 3 , it yields <p 2 = 60° and 4 = 240°. Then, using the value FJ 1 = -0.85 fiC/cm 2 , 
obtained in the DE limit (see Fig. [2]), P cl can be estimated as —0.43 /iC/cm 2 , which is 
consistent reasonably well with P cl = —0.55 /zC/cm 2 , obtained for the noncollinear magnetic 
ground state of YMn03 without additional approximations.- In fact, the main discrepancy 
is caused by the DE limit for PJ. 1 . For example, if one uses PjJ = —1.04 /xC/cm 2 , obtained 
without the DE approximation,- and the angular dependence of P el , given by Eq. ( 1201 . P el 
can be estimated as —0.50 /iC/cm 2 , which is much closer to P cl = —0.55 /iC/cm 2 . 

D. Relative directions of electronic and ionic polarization, and problems of struc- 
tural optimization in LDA+[7 

So far, we considered only electronic polarization, which was induced by the orbital 
ordering in the FM zigzag chains. In this section, we will discuss how this electronic part is 
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related to the ionic polarization in the noncentrosymmetric P1\nm structure. 

Moreover, we will elucidate the microscopic origin of the "order of magnitude difference" , 
which typically exists between experimental and theoretical values of the FE polarization, 
reported for the orthorhombic manganites with twofold periodic magnetic texture. The 
problem is formulated as follows. The great advantage of the first-principles calculations 
is that they allow us to perform the structural optimization and to find theoretically the 
atomic displacements, which are caused by the exchange-striction effects in the E-type 
AFM phase. If one does such structural optimization for the orthorhombic manganites 
and subsequently calculates the FE polarization, the latter will be of the order of several 
/iC/cm 2 .- The conclusion is rather generic and was obtained for several popular types of the 
exchange-correlation functionals, such as LSDA (Ref. s|) , generalized gradient approximation 
(GGA, Refs^ andQ), and LDA(GGA)+£7 (Refs. 5 and[8|). The experimental polarization is 
typically smaller than 0.5 /zC/cm 2 .- On the other hand, if one takes the experimental P2inm 
structure and calculates the FE polarization, it will be at least of the same order of magnitude 
as the experimental one.- 1 ^- The reason of such discrepancy is that, in the experimental 
P2\nm structure, there is a large cancelation of electronic and ionic contributions to the 
FE polarization, while in the theoretically optimized structure, these two contributions have 
the same sign and the cancelation does not occur.- 

In this section, we will further clarify the situation. In orthorhombic manganites, there 
are three types of atomic displacements, which control the FE polarization: 

(i) The Jahn- Teller distortion, which gives rise to the orbital ordering; 

(ii) The exchange striction, which specifies the type of the ordering in the FM zigzag chain 
and, therefore, the sign of the electronic polarization. Note, that in the cenrosymmetric 
Pbnm structure, the FM chains with the 3x 2 —r 2 /3y 2 —r 2 and 3y 2 — r 2 /3x 2 — r 2 type of 
the orbital ordering are equivalent as they build two degenerate magnetic states. This 
degeneracy is lifted in the P2\nm phase by the exchange striction effects, which pick 
up only one type of the FM zigzag chains (characterized by larger Mn-Mn distances). 
As soon as the FM chains are selected, the type of the orbital ordering is fixed, so as 
the sign of the electronic polarization. 

(iii) The FE atomic displacements, which occur in response to the magnetic inversion 
symmetry breaking and control the sign of the ionic polarization. 
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The goal of this section is to understand how these three types of the lattice distortions 
correlate with each other in the experimental and theoretically optimized P2\nm structures 
of YMn0 3 . 

Let us consider the ionic polarization and concentrate on the behavior of the oxygen sites, 
which are located in the ab plane and give the largest contribution to P^ n .— In principles, 
one can consider the contributions of other atomic sites, which do not alter the conclusions. 
Then, P^ n can be presented in the following form: 



where Zj are the atomic charges and At» are the atomic displacements away from the cen- 
trosymmetric positions. Moreover, it is understood that around each Mn site in the primitive 
cell, the summation runs over four oxygen sites, located in the nearest neighborhood of Mn. 
Since each oxygen is shared by two Mn atoms, this leads to the additional prefactor 1/2. 
There are many possibilities for choosing the centrosymmetric reference point for evaluation 
of ATj. The final result should not depend on this choice. For our purposes, it is convenient 
to choose ATj = tq—tm-h (in the other words, we assume that in the centrosymmetric struc- 
ture, all oxygen sites "fall" on the central Mn site). This can be done because Mn sites do 
not contribute to the FE polarization of the ionic type along the orthorhombic a axis.- The 
reason is that, apart from a constant shift, the projections of Mn sites onto the a axis are 
either or a/2 (modulo the lattice translation a) and, therefore, can be transformed to each 
other by the reflection a — > —a. The Mn sites do contribute to the ionic polarization in the 
be plane. However, all these contributions have antiferroelectric character and cancel out 
after summation over the primitive cell. Thus, around each Mn site, the evaluation of P^ n 
is reduced to the summation of At^ over neighboring Mn-0 bonds with the perfactors given 
by Eq. (123]) . Such a construction is very convenient, because in the centrosymmetric Pbnm 
structure, each Mn site is located in the inversion center. Therefore, the sum of At^ over 
all neighboring Mn-0 bonds will be equal to zero. In the Pl^nm structure, however, such a 
construction will give us a finite vector, which can serve as a measure of noncentrosymmetric 
atomic displacements around each Mn site. For our purposes, only the FE (a) components 
of these vectors are important, while the b and c components are antiferroelectric and will 
cancel each other. Using this construction and taking the ionic value Zq = — 2|e|, the con- 
tribution of the planar oxygen sites to P^ n in the experimental P2 1 nm structure can be 




(23) 
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FIG. 9. (Color online) Directions of ionic contributions to the polarization, caused by ferroelectric 
displacements of oxygen atoms around each Mn site in the ab plane of noncentrosymmetric P2\nm 
phase of YMnO"3 (a and c), and the orbital ordering, realized in the ferromagnetic zigzag chain (b 
and d), as obtained for the experimental (a and b) and theoretically optimized structure (c and d). 

estimated as 0.73 yuC/cm 2 , which is totally consistent with the previous finding.- 

The distributions of such vectors, obtained for the experimental and theoretical structures 
of YMnC>3, are shown in Figs. E(a) and[9]^c), respectively. As for the theoretically optimized 
structure, we use results of LDA+£7 calculations with U— 2.2 eV (see Ref. M). Nevertheless, 
we would like to emphasize that very similar results were obtained in LSDA and LDA+£7 
with U— 6.0 eV.- As is seen in Fig. [91 the FE displacements have the same direction in 
the experimental and theoretically optimized P2 x nm structure of YMnC>3. This direction 
corresponds to the positive value of P^ n . 

Corresponding orbital ordering, realized in the FM chains, is shown in Figs.|9]^b) and[9^d), 
for the experimental and theoretical structure, respectively. For the experimental P2inm 
structure, the orbital ordering is of the 3x 2 —r 2 /3y 2 —r 2 type. Therefore, the electronic po- 
larization is negative, and there is a partial cancelation of the electronic and ionic terms, 
which explains a relatively small value of the experimental polarization.- However, the theo- 
retical optimization of the crystal structures, performed both in LSDA and LDA+C/, yields 
different type of the orbital ordering: 3y 2 —r 2 /3x 2 —r 2 instead of 3x 2 —r 2 /3y 2 —r 2 . Therefore, 
the electronic polarization will be positive, and the cancelation does not occur. 

Thus, the directions of FE displacements, obtained in LSDA and LDA+U, are inconsis- 
tent with the type of the orbital ordering, realized in the FM zigzag chains. This seems to be 
a serious problem of the first-principles calculations and at the present stage it is not clear 
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how it should be solved. On the computational side, many attention recently is paid to the 
screened hybrid functionals (see, e.g., Ref. 1281). Therefore, it would be interesting to see how 
these functionals will work for the structural optimization in multiferroic compounds, where 
the inversion symmetry is broken by the magnetic degrees of freedom. The first applications 
for HoMnC>3 seem to show that the problem persists: although the electronic polarization is 
decreased, mainly due to the increase of the on-site level splitting, it has the same sign as the 
ionic one and the total polarization is overestimated in comparison with the experiment.— 
On the other hand, the directions of FE displacements can be controlled by the relativistic 
SO interaction, which is typically ignored in the process of structural optimization. This 



point of view was proposed, for example, in Ref. 



30 



E. Switching electric polarization by changing the magnetic texture 

What is interesting about the multiferroic systems is that the value and the direction 
of the FE polarization depend on the magnetic texture and, by changing this texture, one 
can also change the vector of polarization. In this section, we will discuss how such a 
behavior can be realized in the twofold periodic magnetic texture. Again, let us consider 
the centrosymmetric Pbnm structure and assume that the inversion symmetry is broken 
exclusively by the magnetic order. In such a case, most of attention is focused on the E- 
type AFM phase (Fig. EJ), which breaks the inversion symmetry but preserves the symmetry 
operation {C^\a/2+b/2}. Therefore, the FE polarization will be parallel to the a axis. 

Now, the question is whether there are other types of the magnetic texture, which would 
break the inversion symmetry. As an example, let us consider the magnetic texture in 
Fig. HPT a). In the plane z = 0, it is identical to the £7- type AFM order, and can be trans- 
formed to itself by applying the symmetry operation {C^\a/2+b/2} around even magnetic 
sites 2 and 4. Alternatively, one can apply the symmetry operation {C*;^— a/2— 6/2} around 
odd magnetic sites 1 and 3. In the .E-phase, the same symmetry operations can be applied 
in the planes z = ±c/2 and also will transform the plane z = c/2 to the equivalent to it 
plane z = —c/2. The magnetic texture in Fig. [TuTa) is obtained by the additional inversion 
around odd magnetic sites in the plane z = c/2, which interchanges the symmetry opera- 
tions {C%\a/2+b/2} and {Cl\— a/2— b/2}. Thus, the plane z = c/2 can be transformed to 
itself by the symmetry operation {C^\a/2+b/2} around odd sites or by {C%\—a/2—b/2} 
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<*Kdeg) 

FIG. 10. (Color online) (a) Antiferromagnetic texture yielding finite ferroelectric polarization along 
the c axis, (b) Behavior of electronic polarization in YMn03 upon rotation of magnetic moments, 
as obtained in the double exchange model with the Hartree-Fock potential Vj (DE +U). In the 
rotated texture, the directions of spins at the sites 1 and 3 were fixed, while the spins at the sites 
2 and 4 were rotated by the angle 4>, as explained in panel (c). The interlayer coupling was kept 
AFM for the sites 1 and 3 and FM for the sites 2 and 4. 

around even sites. Therefore, the symmetry operations {C%\a/2+b/2} and {C^|— a/2— 6/2}, 
although preserved locally in each of the plane, are broken globally, because they cannot 
simultaneously transform the planes z = and z = ±c/2 to themselves. 

Instead, the magnetic texture in Fig. [TUf a) obeys the symmetry operation {C^\c/2}, 
which is another symmetry operation of the space group Pbnm. Therefore, the FE polar- 
ization in this phase will be parallel to the c axis. According to the above arguments, each 
plane may carry a finite polarization parallel to the a axis. However, since neighboring 
planes are connected by the symmetry operation {C^\c/2}, the contributions from different 
planes will cancel each other. 

The behavior of P||c, obtained in the DE model for YMnC>3, is explained in Fig. [TOT b). 
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Here, we again consider a continuous rotation of spins between two kinds of the AFM 
domains via an intermediate spin-spiral phase, as explained in Fig. ITOT c). In comparison 
with Fig. [2J the planes z = and z = c/2 are connected by the FM pathes between even 
magnetic sites. 

P||c appears to be about two orders of magnitude weaker than P||a in the .E-phase 
(Fig. [2]). Nevertheless, this result is very natural and can be easily understood by considering 
the perturbation theory arguments, similar to the ones in Sec. Ill I CI Namely, in order to 
obtain P| |c, we should consider the transfer integrals ty between all possible combinations of 
sites i and j along the c axis. Of course, the main contribution is expected from the NN sites. 
Moreover, according to Eq. ( 122]) . in order to contribute to P| |c, these transfer integrals should 
have both symmetric and antisymmetric components. However, due to the combination of 
{C 2 \c/2} and the inversion symmetry around the Mn sites, the NN integrals between the 
planes z = and z = c/2 will satisfy the following property: R 2 c iij{R 2 c ) T = tji, where the 
matrix transformation R 2 , corresponding to the 180° rotation around the c axis, changes 
the sign of some of the matrix elements of tij. Therefore, in the crystal- field representation, 
one can always choose the phases of the basis orbitals such that the corresponding matrix 
of the transfer integrals ty would become totally symmetric. Thus, the NN contributions to 
P||c in the second order of 1/A will vanish, and P||c has finite value due to either next-NN 
integrals, which are small (all transfer integrals for YMn03 can be found in the Supplemental 
Material of Ref. |8[) or the higher-order effects with respect to 1/A, which are also small. 
This naturally explains the fact that P||c is much smaller than P||a. 

This finding resembles the behavior of multiferroic manganites with nearly fourfold pe- 
riodic magnetic texture, for which the possibility of switching the electric polarization was 
demonstrated experimentally.-^ For example, in TbMnOs the polarization is aligned along 
the orthorhombic c axis. However, the external magnetic field applied along the b axis 
will change the magnetic texture and align the polarization parallel to the a axis.- More- 
over, most of experimental data also confirm the fact that P||c is smaller than P||a. For 
example, such a behavior is typical for the Eui-^Y^MnOs compounds, containing only non- 
magnetic rare-earth elements, that excludes the influence of the 4/ magnetism on the FE 
polarization.-^ The results of this section suggest that this behavior is more generic and 
can be anticipated in other regimes, including the twofold periodic magnetic systems. The 
origin of this phenomenon is related to the specific symmetry of the crystal structure (in 
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the case of orthorhombic manganites - the Pbnm symmetry) and how it is lowered by the 
magnetic ordering in the DE limit. It should not be confused with the spin-spiral alignment, 
which does break the inversion symmetry of the DE Hamiltonian (see Sec. IIII CD . 

In is interesting to note that the magnetic texture depicted in Fig. [TO] can be viewed 
as a "defected E-type AFM texture", where the "defects" are two FM bonds between the 
planes z = and z = c/2. Of course, such "defects" are energetically unfavorable and, after 
including the SO interaction, this magnetic texture will change in order to minimize the FM 
coupling in the defected bonds. This will lead to the substantial deformation of the mag- 
netic texture in Fig. HUT a). Nevertheless, we would like to emphasize that the noncollinear 
magnetic texture with P||c can be stabilized even after including the SO interaction. The 
situation was discussed in Ref. LZJ- 



IV. DISCUSSIONS AND CONCLUSIONS 

This work is a continuation of previous studies, devoted to multiferroic manganites, which 
crystallize in the orthorhombic Pbnm and P2inm structure.- 1 ^ Our main motivation was to 
present a transparent physical picture, which would explain why and how the ferroelectric 
polarization is induced by some complex magnetic order. For these purposes we invoke 
the double exchange theory, which was formulated for the low-energy model, derived from 
the first-principles electronic structure calculations. As far as the polarization is concerned, 
the DE theory is very robust and reproduces results of more general mean-field Hartree- 
Fock calculations at a good quantitative level. Furthermore, the main advantage of the DE 
theory is that it allows us to greatly simplify the problem and, in a number of cases, derive 
an analytical expression for the FE polarization. Thus, we could clarify very basic aspects 
of the FE activity in manganites with twofold periodic magnetic texture. 

In our analysis we started from the general Berry-phase theory.— 1 ^ In the case of im- 
proper ferroelectrics, the basic quantity to be considered is the electronic polarization, which 
incorporates the change of the electronic structure in response to the noncentrosymmetric 
alignment of spins. Then, our main message is that, for the analysis of electronic polar- 
ization in realistic manganites, one can always use two physical limits. The first one is 
the limit of large intra-atomic splitting A ex between the majority- and minority-spin states. 
The second one is the limit of large intra-atomic splitting A between the majority-spin e g 
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states. Therefore, for the electronic polarization, one can always consider the perturbation 
theory expansion with respect to both 1/A ex and 1/A. This perturbation theory describes 
asymmetric transfer of some weight of the Wannier functions to the neighboring sites, which 
gives rise to the polarization. 

There is some similarity with the theory of superexchange interactions, which deals with 
the virtual hoppings,™ and where the terms proportional to 1/A and 1/A CX account for 
the FM and AFM contributions, respectively.— Therefore, the DE limit A ex — t-oo would 
correspond to neglecting all AFM contributions. It may not be a good approximation 
for interatomic magnetic interactions. Nevertheless, the main difference for the electronic 
polarization is that it appears only in the second order with respect to 1/A and 1/A ex . The 
physically relevant picture corresponds to the situation where A ex > A. Then, due to the 
inequality (A/A ex ) 2 <C A/A ex , it is logical to keep the effects of the first order of 1/A ex 
in the analysis of superexchange interactions, but neglect the effects of the second order of 
l/A^ i n the analysis of electronic polarization. This again justifies the use of the DE limit 
in the latter case. 

On the basis of this perturbation theory expansion, we were able to explain how the 
electronic polarization depends on the relative directions of spins in an arbitrary twofold 
periodic magnetic texture. Particularly, the multiferroic effect in orthorhombic manganites 
is a nonlocal phenomenon in the sense that the inversion symmetry is broken by making 
some of the Mn-Mn bonds magnetically inequivalent. In the DE model, this inequivalence 
is achieved by the additional modulation of transfer integrals by Then, one trivial 
conclusion is that there will be no magnetic inversion symmetry breaking in the spin-spiral 
phase, where all £y are the same. Therefore, in order to make finite polarization, it is essential 
to deform the spin spiral. In orthorhombic manganites, such deformation is caused by the 
relativistic spin-orbit interaction.-^ The second important precondition for the FE activity 
is the asymmetry of the transfer integrals, which should simultaneously have symmetric and 
antisymmetric components. 

We also pointed out on a serious problem in the structural optimization, which appar- 
ently exists in the first-principles calculations (at least at the level of LDA+£7 and GGA+U 
approximations for the exchange- correlation functional without relativistic spin-orbit cou- 
pling) and which typically results in the large overestimation of the value of FE polarization 
in comparison with experimental data.- In this work, we were able to clarify the origin 
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of this problem: in the theoretical structure, the directions of noncentrosymmetric atomic 
displacements are inconsistent with the type of the orbital ordering in the ferromagnetic 
zigzag chains, which controls the sign of the electronic polarization. As the result, the elec- 
tronic and ionic contributions have the same sign in the theoretically optimized structure, 
while, according to the experimental crystal structure, they should have opposite signs and 
partially cancel each other. 

Finally, we explained how the electronic polarization can be switched between orthorhom- 
bic a and c directions by inverting the magnetic texture in every second ab plane. We also 
expect a gigantic change of the absolute value of the polarization itself, which is related 
to very different symmetry properties of the nearest-neighbor transfer integrals along the c 
direction and in the ab plane of manganites. 

In this work, our analysis was limited by twofold periodic magnetic textures, which 
illustrate the basic idea of the double exchange theory of ferroelectric polarization. The idea 
can be extended to the systems with more general magnetic periodicity: apart from the 
additional complexity of the magnetic texture, there is no fundamental difference between 
twofold and more general magnetic periodicity. In both cases, the basic property, which 
should be considered and which gives rise to the ferroelectric activity is the alternation of 
angles between spins in different Mn-Mn bonds. 
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